function f = zifilter(order, bw3dB)
%ZIFILTER Summary of this function goes here
%   Detailed explanation goes here
    s=tf('s');
    tau_N = ziBW2TC(bw3dB, order);
    f = 1/(1+tau_N*s)^order;
end

